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Abstract 

Kinetic energy spectra and fragment velocity correlations, simulated by means of 
stochastic mean-field calculations, are successfully confronted with experimental 
data for single multifragmenting sources prepared at the same excitation energy 
per nucleon in 32 AMeV i29xe+"'^*Sn and 36 AMeV Gd+'^''^\] central colli- 
sions. Relying thus on simulations, average freeze-out times of 200-240 fm/c are 
estimated The corresponding spatial distributions of fragments are more compact 
for the lighter system (~3-4Vo vs ~8Vo). 
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1 Introduction 



A detailed investigation of multifragmentation in heavy-ion collisions at inter- 
mediate energies is of great interest at present time in connection with phase 
transition in finite nuclear systems [1,2,3,4,5,6,7]. Fused systems produced in 
central collisions between heavy nuclei allow to address rather large pieces 
of nuclear matter, of about 200-400 nucleons which undergo multifragmen- 
tation [8,9]. The present paper, presenting intra-event correlations which are 
highly enlightening for the origin and the features of the process, enlarges and 
completes the comparison started in reference [9]. 

Two nuclear reactions, at similar available energy per nucleon were exper- 
imentally studied: 32 AMeV i29Xe+"«*Sn and 36 AMeV i55Gd+"»*U. Single 
multifragmenting sources at the same excitation energy per nucleon are formed 
in central collisions. Their properties were examined in detail [9,10,11,12]. The 
angular and average energy distributions of all fragments {Z > 5), isotropic in 
the cm., are compatible with those of a thermalised multifragmenting source. 
The two measured average fragment multiplicities are in the ratio of the to- 
tal charges of the systems, while the charge distributions are identical [13]. 
This experimental observation can be taken gnature either of a large 

exploration of phase space or of volume instabilities. Up to now all global 
reaction characteristics, except the widths of the individual or total kinetic 
energy distributions of charged products, were equally well described by sta- 
tistical [10,14,15,16] and dynamical [9] simulations. The two interpretations 
are not contradictory since one can suppose that the dynamics of the collisions 
is sufficiently chaotic to explore enough of the phase space, allowing a statis- 
tical description of the fragment production. Higher order charge correlations 
generated in the same dynamical approach were also successfully confronted 
with the data [17,12] for the 32 AMeV ^29Xe+"»*Sn case. Both calculated 
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and experimental results suggested an enhanced production of equal-size frag- 
ments, compatible with a multifragmcntation induced by homogeneous density 
fluctuations in the liquid-gas coexistence region (spinodal decomposition). 

We continue a step further by examining detailed properties of fragment ki- 
netic energy spectra and their intra-event correlations. The validity of the 
dynamical simulation results versus the experiment will be tested through 
these observables. Fragment velocity correlations are supposed to bring space 
and time information concerning the multifragmentation source [18] and may 
thus put supplementary constraints on models. They may allow to trace back 
the event topology at "freeze-out" , when the fragments become free and feel 
only the Coulomb repulsion. 

The method and the data are described after the experimental selection of 
events and the detailed presentation of kinetic energy properties of fragments. 
Then results of dynamical simulations are compared to the experimental pat- 
terns. The consequences on the fragment spatial distribution at the end of the 
elapsed time corresponding to their separation are discussed too. Conclusions 
are finally drawn. 



2 Experimental selection and kinetic energy spectra 

The data were collected, as largely described in ref. [8], by means of the 336 
multilayer detection cells of the An multidetector INDRA [19]. Accurate frag- 
ment identification and energy calibration - crucial for this type of studies - 
were achieved by taking into account pulse height defects in the silicon de- 
tectors [20] and the influence of the quenching and of the 5-rays on the light 
output of the Csl(Tl) scintillators [21,22]. The energy of the detected prod- 
ucts is obtained with an accuracy of 4%. In the following, results relative to 
the kinetic properties of fragments are discussed in the centre of mass of the 
reactions and thus depend on the masses (not measured) attributed to the 
fragments. In all figures, a single mass close to that of the /3-stability valley 
was attributed to each fragment of charge Z. However recent studies on the 
de-excitation of hot heavy nuclei have shown that the cold residues have an 
average mass lying on the evaporation attractor line [23], corresponding to 
neutron deficient nuclei. The two mass formulae become different for nuclei 
with charges larger than 20-25; thus the chosen relation is only important for 
the heaviest fragments discussed in the present paper. It was verified that the 
average cm. energy of these fragments changes at most by a few MeV when 
changing the hypothesis on their mass, which do not alter the conclusions 
given below. 

Among the complete experimental events: total pseudo-momentum >75% 
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of the entrance channel value {ZprajV proj) , and total detected charge >80 
(Xe+Sn) or >120 (Gd+U), compact single sources were selected with the 
constraint of flow angle > 60° [24,25,8,9]; indeed it was shown in previous 
studies that while events present the topology of emission from two sources 
at small flow angles, they evolve towards a single-source configuration above 
60° (see figure 9 in ref [8] and figure 1 in ref [6]). A minimum number of 
three fragments (Z>5) was additionally required. The measured average frag- 
ment multiplicities are 4.3 for Xe-|-Sn and 6.2 for Gd+U [9]. The reconstructed 
sources of both systems have an excitation energy per nucleon ~ 6.5 MeV, and 
average charge and mass (estimated following the method described in [9]): 
~99, ~236 and ~150, ~378: the missing charge, mass and energy relative to 
the entrance channel values are carried away by high energy particles (Z<2) 
emitted backward and forward. 

The cm. kinetic energy spectra of various fragments emitted in these central 
collisions are shown in fig. 1. They present asymmetric shapes for the lighter 
fragments, and tend to become more symmetric for Z>15, encoding informa- 
tion about the Coulomb repulsion, radial expansion and the temperature of 
the single source from which they originate. 

More information on the kinematical characteristics of the fragments was ac- 
cessed by sorting events according to the fragment multiplicity, M/, and the 
rank of the fragment in the event (largest, Zmax, second largest, Zmax2, third 
largest, Z^axs, and so on). An example of the results is shown, for the Xe+Sn 
system, by the symbols in fig. 2. Events with four fragments are displayed, 
but we have observed that the experimental patterns and energy values are 
identical irrespective of the fragment multiplicity. Ref. [26] shows that, around 
30AMeV, the fragment kinetic energy is essentially from Coulomb origin; this 
invariance is a good indication that events arise from sources with very similar 
charge, independently of the fragment multiplicity. The average kinetic energy 
of fragments first increases with the fragment charge, and then saturates and 
even decreases for charges Z>20-25, independently of the fragment multiplic- 
ity. This evolution is a Coulomb effect. The most striking feature in fig. 2 is the 
particular role played by the largest fragment in each partition: in the region 
Z— 15-25 the heaviest fragment, Z^nax, has always the lowest average kinetic 
energy. Note that this behaviour was already observed for the same system 
at 50 AMeV in ref. [25]. This effect may be caused by an inhomogeneity of 
the system created by cm. conservation constraints, positioning the largest 
fragment close to the centre [27]; thus Coulomb as well as radial expansion 
energies would be reduced for this largest fragment. 

The same picture was drawn for the heavier system, Gd+U. Again one ob- 
serves the same pattern whatever the fragment multiplicity. Because of the 
smaller number of events collected for this system, fig. 3 displays results with- 
out multiplicity selection. As expected from Coulomb effects, the average en- 
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Fig. 1. Kinetic energy spectra of a few fragments for 32 AMeV Xe+Sn - top panel 
- and 36 AMeV Gd+U - bottom panel. Dynamical simulation results (histograms) 
are compared to experimental data (symbols). 



ergies are larger for Gd+U than for Xe+Sn. Here again the average energy of 
the largest fragment is smaller than the energy of the other fragments, and 
shows a maximum for Z~30-35 [9]. However the relative difference between the 
energies of the largest fragment and those of the others is smaller for Gd+U 
than for Xe+Sn. 
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Fig. 2. Average cm. kinetic energy of fragments in events with Mj=4, for the Xe+Sn 
system at 32 AMeV. Symbols show the experimental values and the lines the results 
from the BOB simulation. A minimum of 10 fragments of given Z was required for 
each plotted point. 
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Fig. 3. Same as figure 2 for the Gd+U system at 36 AMeV and all fragment multi- 
plicities. 

The evolution of these average energies gives first information on the topology 
of the events when nuclear interaction becomes ineffective; a connected ob- 
servable is the average fragment charge as a function of their multiplicity and 
rank, as shown in fig. 4. Only Z^ax strongly varies with Mj; the fragments of 
rank larger than 2 keep almost constant charge values. The largest fragment 
again presents a special behaviour. The observation of figs. 2, 3, 4 reminds 
of the classification of multifragmentation events [28] in "sun" (a large cen- 
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Fig. 4. Measured average charges of the five largest fragments of events versus the 
fragment multipUcity. The vertical line represents the average fragment multiplicity. 

tral fragment surrounded by small ones) for small Mj- events and "soup" (a 
random spatial distribution of fragments of comparable sizes) for large M/ 
events. However, in this last case, one would expect that the largest fragment 
looses its specificity, which is not true here as it keeps a smaller average en- 
ergy, even at large multiplicities. This above image is thus too simple, and 
one may hope to gain a deeper understanding from the velocity correlations 
presented in the next section. Several advanced analyses took advantage of 
the special behaviour of the largest fragment, for instance universal fluctua- 
tions [29] which indicate a transition from an ordered to a disordered phase, 
and bimodality [30,31] which suggests the occurrence of a first order phase 
transition in the studied systems. 



3 Velocity correlations 

Relative velocity correlation functions between emitted products are tools to 
get information on emission time scales and on the disassembling source vol- 
ume. Light charged particle correlation functions were heavily investigated 
in the past (see [32,33] for reviews and [34,35,11]). A new technique, tak- 
ing advantage of both the height and the width of the correlation functions 
was recently developed, allowing to disentangle fast and slow emission com- 
ponents [36]. Fragment correlation functions were obtained for very light frag- 
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ments [37,38], and over a much larger charge range [39,40]. In this case energy 
and momentum conservation becomes important. The Coulomb interaction 
between the two fragments, and between the fragments and the source starts 
to dominate; one may thus expect to get information on the arrangement 
of the fragments inside the source volume [41]. Fragment velocity correlation 
functions are often presented as a function of a reduced velocity, Vj-cd, in order 
to increase the statistics, through mixing of fragments with different charges 
by scaling the Coulomb effects. 

Fragment velocity values are derived from their energy using estimated masses 
(see sect 2). The emission polar and azimutal angles of a fragment are chosen 
randomly over those covered by the module which detected this fragment. For 
each couple i, j of fragments with velocities Vi and i*], in an event of multiplicity 
M/, (i, j = 1, Mf and i ^ j), the reduced relative velocity is defined as [18]: 

^^red OC , = , (1) 



The formalism for two particle correlation depending on one variable was 
considered for intermediate mass fragments (IMF) emitted in multifragmen- 
tation [41,42,28]. We define the two fragment correlation function dependent 
on v^ed as: 

^{vM{vuZ^K^,=const 111 [(^^k, ^k)] [(^^1, Zl)] 



where n2[(^^! ^i)? (^^! ^j)] is the probability to find two fragments of atomic 
numbers Zi, Zj with the reduced velocity v^ed in one event, while the product 
ni[('J\, -^k)] ni[('j^) Zi)] is the probability to find two fragments with the same 
reduced velocity, but in two different events. The sum at the numerator in 
eq. (2) spans all couples having v^-cd = const in all real events and it accounts 
for correlated fragments, while the denominator accounts for the uncorrelated 
ones. The reduced velocity distribution function of the correlated fragments 
is built by taking into account all two-fragment combinations: C{Mf,2) — 
Mf\/{2\{Mf — 2)!) in one event and all the experimental selected events. For 
the uncorrelated case, we have proceeded as follows. For a given real event, 
twenty pseudo-events were generated by associating to each fragment another 
one with the same charge, randomly found in one of the other experimental 
events of the same sample. The fragment multiplicity and the sum of the 
fragment charges of the initial event are thus conserved. A pseudo-event is 
validated only if all its fragments come from different real events and were 
detected in different INDRA modules. This procedure reduces the biases due to 
Coulomb effects. The ratio of the two distribution functions: real over pseudo- 
events, in each reduced velocity bin, leads to the correlation function defined in 
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eq. (2). No normalisation factor but the factor 20 coming from the number of 
uncorrelated events is used. Since only fragments from the selected sample are 
used to construct the uncorrelated events, their reduced velocity distribution 
essentially reflects the two-fragment phase-space population of the detection 
array in the absence of any final state interaction. 

Several correlation functions are proposed here, depending on the size of the 
fragments considered in the procedure: 

i) all fragments considered (Zj j > 5); 

ii) intermediate mass fragments (IMF) 5 < Zij < 20; 

iii) reduced velocity correlation of the heaviest fragment Z^ax with each of the 
others in the event Z\ < Z^ax- 

In the absence of correlation the value of the correlation function should be 
equal to 1. At small reduced velocities momentum conservation as well as 
Coulomb repulsion influence the correlation functions. 

The three types of experimental correlation functions are shown as symbols 
in fig. 5, for 32 AMeV Xe+Sn on the left side and for 36 AMeV Gd+U on the 
right side. It was verified with the help of the simulation described in the next 
section that the correlation functions are negligibly affected by the granularity 
of the INDRA array. In the six panels one observes a dip in the vicinity of 
■^^red = 0, followed by an enhancement in the intermediate reduced velocity 
region (0.015-0. 035c). The width of the hole, independent of the charges of 
the two fragments considered (see fig. 5. 14 of ref. [10]), refiects the distances 
between fragments at freeze-out. Its depth, with a minimum value very close 
to zero, indicates that fragments are not emitted independently, as assumed 
in sequential decay models [28]. Indeed it was shown in [43], for a system 
similar to Xe+Sn, that the Coulomb interaction effect for two fragments is 
negligible when the time between their emission is long: ~ 10"^ fm/c. The 
height of the bump at intermediate t'red is expected to be connected to the 
distribution of the nuclear matter inside the source volume at freeze-out, for 
instance, as mentioned above, "sun events" - evidenced by a very pronounced 
peak, or "soup events" for which the peak is diminished and the distribution 
is fiatter [28]. 

The shape of the correlation function evolves, for each system, with the size 
of the correlated fragments. The height of the bump at intermediate v^cd is 
small when all the fragments are considered - case i), it becomes fiatter when 
only small fragments are treated - case ii). Then, it increases - case iii) - when 
the heaviest fragments are isolated in the procedure. 

In each of the three cases, the width of the hole at a value of the correlation 
function equal to 0.5 is very similar for the two systems. The shortest intra- 
fragment distances seem thus to be quite independent of the system size. For a 
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Fig. 5. Fragment-fragment correlation functions for 32 AMeV Xe+Sn - left panels - 
and 36 AMeV Gd+U - right panels. Dynamical simulations (histograms) are com- 
pared to experimental data (solid symbols). When not visible, the statistical error 
bars are smaller than the size of the symbols. Errors on the simulated functions are 
similar to the experimental errors for the Xe-|-Sn function of the same type. See 
text for explanations of cases i)-iii). 

given system, the width of the Coulomb hole is slightly narrower and the peak 
more pronounced for cases iii), testifying about a stronger Coulomb interaction 
of the heavy fragments with the others, in relation with their position in the 
source. This effect is especially visible for the system Xe+Sn. In fact, a more 
pronounced peak for this system is manifest in all three cases, proving may 
be simply that, at smaller average multiplicities, the centre of mass is closer 
to the heaviest fragment and the distances between this fragment and all the 
others are on average shorter. While for cases i) and ii) it was verified that 
the correlation function remains at 1 up to the maximum range populated by 
the reduced velocity distribution (~0.09c), for case iii) the correlation func- 
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tion strongly decreases below 1 above ~0.04c, particularly for Xe+Sn. Indeed 
the range of the reduced velocity distributions (both correlated and uncorre- 
lated) are shorter in that case (~0.05c); it was claimed in [42] that energy 
and momentum conservation implies that at very large reduced velocities, the 
correlation function will tend to zero; it is not clear why this trend would ap- 
ply only to correlation functions of type iii). The reason of the decrease must 
probably be found in the uncorrelated distributions; it is indicated in [40] 
that the shape of the correlation function is very sensitive to the background 
constraints when the heavier fragment is isolated. 

Quantitative information on the source sizes and fragment emission times will 
be derived from the confrontation with full dynamical simulations of collisions 
in the next section. 



4 Compcirisons with collision simulations 

4-1 BOB simulations 

One way to explain multifragmentation in heavy ion collisions at Fermi en- 
ergies is to correlate it with a phase transition of liquid-gas type, due to the 
specific form of the nucleon-nucleon interaction. After a compression stage, 
the nuclear system enters an expansion phase, during which it cools down 
and evolves in the spinodal region (of negative compressibility) of the phase 
coexistence domain, where multifragmentation occurs through the growth of 
density fluctuations [44]. Stochastic mean field simulations of nucleus-nucleus 
collisions, based on the Boltzmann-Langevin equation, allow for the treat- 
ment of unstable systems [45,46,47]. Nevertheless applications to 3D nuclear 
collisions are prohibited by computational limitations. The dynamical path 
through the spinodal region has been instead simulated by means of a Brow- 
nian force [48,49] - Brownian One Body (BOB) dynamics - grafted, at the in- 
stant of maximum compression, ~ 40 fm/c, on the one-body density evolution 
calculated in a Boltzmann-Nordheim-Vlasov (BNV) approach. The chosen 
self-consistent mean field potential gives a soft equation of state. The ingre- 
dients of the BOB simulations are presented in detail in [9], as well as the 
comparison between average observables concerning filtered simulated events 
and experimental ones. In the simulations of head-on collisions, both systems 
form a single source which expands with time and breaks-up into several frag- 
ments (as shown by the time evolution of the densities, figure 1 of ref [50]). 
At 100 fm/c, at low density inside the spinodal zone, the systems are already 
thermalised, with a temperature of 4 MeV and a small radial velocity at the 
surface (~ 0.1c), and the first fragments appear. An algorithm for reconstruct- 
ing fragments is applied at intervals of 20 fm/c, based on minimum density 
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Fig. 6. a) Evolution with time of the average fragment [Z > 5) multiphcity calcu- 
lated for Xe+Sn and Gd+U collisions [9]; b) variation rate of the same quantity. 

for each system, when the average fragment multiphcity becomes a constant, 
as shown in fig. 6a). For a given system, this multiplicity is independent of 
the value of pmin [9] • The variation rate of the fragment multiplicity shows in 
fig. 6b) that any dramatic change in the average number of fragments ends at 
~ 200 fm/c for ^29Xe+""*Sn and at ~ 240 fm/c for i^^Gd+^^^U. These values 
can be considered as reasonable average freeze-out times [50]. 

At the end of a BOB calculation, the fragments, still hot (e* ~ 3 AMeV), de- 
excite via secondary particle emission, treated by means of the suitable part of 
the code SIMON [51]. Global observables, like Mf, Z, Zmax and Zbound = 
distributions are very well reproduced by BOB calculations, while the average 
fragment kinetic energy is underestimated by about 20% over all the Z range 
for Xe+Sn, and for Z>20 for Gd+U [9]. Part of the discrepancy might come 
from the semi-classical treatment of BOB which underestimates the radial 
energy - by about 50% - as compared to a quantal treatment [7]. 



4-2 Average energies and correlation functions 



After filtering the calculated events through a software replica of the INDRA 
array, the average calculated kinetic energies of fragments were sorted follow- 
ing the fragment multiplicity and rank, as done for experimental data. The 
results are displayed by the lines in figs. 2 and 3. As mentioned above, the 
calculated energies underestimate the experimental values for Xe+Sn while 
they agree with data for the lighter elements in Gd+U. Besides a drawback in 
the calculation mentioned above, this may sign a more compact shape of the 
experimental configuration with respect to the calculated one for Xe+Sn. In 
agreement with experiment the average energy for a given charge is indepen- 
dent of the fragment multiplicity (not shown) and the average energy of the 
largest fragment decreases when its charge becomes larger than ~25 (Xe+Sn) 
or ~30-35 (Gd+U). For the heavier system, the calculated energies show the 
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same hierarchy with the fragment rank, namely the energy of a given charge 
is smaller when the fragment is the heavier of the partition. Conversely, the 
calculated energies are independent of the fragment rank for Xe+Sn. This may 
indicate that the topology of the simulated events reflects well the real one for 
Gd+U while it is slightly different for Xe+Sn. A different topology also influ- 
ences the calculated values of the fragment energy: as stated in [7] Coulomb 
acceleration is more effective for a uniform distribution than for a hollow con- 
figuration as that obtained in BOB [50] . Thus both the larger absolute values 
of the fragment energies, and the smaller energy of the largest fragment may 
indicate that, for Xe-|-Sn, the experimental freeze-out configuration is more 
compact, more uniformly filled, with the largest fragment closer to the centre, 
than the calculated configuration. 

We have also confronted the simulated kinetic energy spectra of the final frag- 
ments with the experimental ones. Fig. 1 shows that, for Xe-|-Sn, the calcu- 
lated spectra are narrower than the experimental ones, and more so for higher 
charges. For the heavier system the agreement between calculated and mea- 
sured spectra is better: both the average values and the widths of the spectra 
are reproduced to within a few percents for odd charges and about 10% for 
even charges, while the differences on these quantities lie between 10 and 20% 
for Xe+Sn. A calculated odd-even effect can be noted: not only is the produc- 
tion of even charge fragment favoured, but both the average values and the 
widths of their spectra are farther from the experimental values than those of 
the odd-charge fragments. Note that for both systems fragments with charge 
Z=5 are strongly underestimated in the calculation. While the odd-even effect 
has to be attributed to the de-excitation process, as no such effect is seen in 
the primary distributions, the disagreement between calculation and experi- 
ment, besides the above mentioned differences in topology, may originate from 

i) the dropping-out of the light charged particles in the last step of the cal- 
culation: only fragments are input to the SIMON code, the particles already 
free (27% of the initial charge for Xe+Sn, 22% for Gd+U) do not participate 
to the Coulomb propagation. This contributes to the underestimation of the 
energy. 

ii) the spatial distribution of these primary particles: placed mainly in the cen- 
tre of the source instead of being uniformly distributed, they would increase 
the Coulomb effect. 

To summarize, these dynamical calculations fit rather well the individual ex- 
perimental energy spectra - a difficult task, scarcely reported up to now for 
other models. 

The method of reduced velocity correlations, as described in the previous sec- 
tion, is applied in the following to fragments from BOB simulated events. A 
good agreement between the simulated (histograms) and experimental (sym- 
bols) correlation functions is observed in fig. 5. The depth and the width of 
the Coulomb hole are well reproduced by the present dynamical simulations 
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for the different correlation functions. Note that the decrease below 1 of the 
correlation functions of type iii) is present also in the simulations. For Gd+U, 
the agreement between the experimental and calculated correlation functions 
is excellent; the small filling- up at low reduced velocity is introduced by the 
filtering step. For Xe+Sn, the three calculated functions are slightly narrower 
than the experimental ones and the peak is displaced towards smaller reduced 
velocities. This discrepancy confirms the difference between the experimen- 
tal and calculated event topology already mentioned. However it should be 
stressed that, for both systems, the correlation functions generated by the 
BOB simulations account better for the experimental ones than statistical 
models like SMM or MMM [10,52,16]. 

4.3 Freeze-out configuration: volume estimate and topology 

The agreement between calculated and experimental kinetic properties, es- 
pecially for the heavier system, encouraged us to give signification to the 
positions of the fragments just after their formation. We analysed the calcu- 




Fig. 7. Fragment positions as a function of the atomic number for: 32 AMeV 
^^^Xe+^^^Sn at 200 fm/c from the beginning of BOB calculation - upper panel, 
240 fm/c - lower panel and 36 AMeV ^^^Gd+'^^^U at 240 fm/c - middle panel. The 
contour scale is logarithmic. 

lated events at the moment when their multiplicity variation rate vanishes: at 
200 fm/c from the beginning of the calculation for Xe+Sn and at 240 fm/c 
for Gd+U (see fig. 6b). At the moment of the fragment separation, the values 
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of their distances ri from the source cm. or of their relative distances In — r^l 
give information concerning the topology of events and associated freeze-out 
volumes. 
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Fig. 8. Radial charge distributions. Upper panel: 32 AMeV ^^^Xe+^^^Sn at 200 fm/c 
from the beginning of the BOB calculation; middle panel: 36 AMeV ^^^Gd+^^^U at 
240 fm/c; lower panel: as in the upper panel but at 240 fm/c. The arrows indicate the 
radii used to estimate the volumes reported in Table 1 (see text for explanations). 

The position radii r of the fragments relative to the cm - fig. 7 - are shorter 
in the Xe+Sn reaction at 200 fm/c (upper panel) than in Gd+U at 240 fm/c 
(middle panel), which may be attributed to the shorter propagation time on 
one hand, and to the Coulomb interaction (Ecoui) and the radial expansion 
i^rad), both inferior for the smallest source [9]. The effect of the time difference 
is shown on the lower panel of the figure. The ridge lines of the plots for the 
two systems have the same slope [53], r decreases when Z increases. These 
observations explain the behaviour of the fragment-fragment correlation func- 
tions shown in fig. 5: the heavier fragments take more central positions, at 
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least at the end of the fragment separation time. In case of Xe+Sn system 
the heavier fragments are moreover very close to the centre of mass (2-3 fm), 
while they are farther away for Gd+U. Correlatively the peak in the correla- 
tion functions increases when the heavier fragments are involved and is more 
pronounced for the lighter system. 

Table 1 

Some characteristics of the 32 AMeV Xe+Sn and 36 AMeV Gd+U single sources - 
calculated by means of BOB dynamical model [9] - in the spinodal region (100 fm/c), 
columns 2 and 3, and at the moment when the mean multiplicity of fragments, Mf , 
saturates. In the last column, the ratio of the average volume of a sphere englobing 
nearly all the fragment centres and the source volume at normal nuclear density. 



system 


Atot 


Ztot 


t (fm/c) 






i29Xe+"9Sn 


238 


100 


200 


5.1 


2.8 


155Gd+238u 


360 


142 


240 


8.1 


7.9 


i29Xe+"9Sn 


238 


100 


240 


5.1 


4.2 
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Fig. 9. Fragment relative distances vs the product of the atomic numbers of the frag- 
ments considered in each couple, for: 32 AMeV ^^^Xe+^^^Sn at 200 fm/c from the 
beginning of BOB calculation - upper panel, 240 fm/c - lower panel, and 36 AMeV 
^^^Gd+^^^U at 240 fm/c - middle panel. The lines show the minimum distances 
between two neighbour fragments. The contour scale is logarithmic. 

The local charge concentration per unit volume, dZ/dV, plotted in fig. 8 as 
a function of the distance r from the cm. for both systems, summarizes well 
what was discussed all along the paper: the configuration for Xe+Sn is more 
compact (and one derives from the experiment a configuration even more 
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uniformly filled) than the one observed for Gd+U. This is also well illustrated 
by the relative distances between two fragment centres Vij = |r^ — rjl, shown in 
fig. 9. The shortest distances concern the closest neighbours, the longest ones 
are representative of the size of the sources. Relative distances are shorter 
for the reaction Xe+Sn at 200 fm/c and also at 240 fm/c - fig. 9 (upper 
and lower panels) than for Gd+U at 240 fm/c (middle panel). The profile 
of each distribution is horizontal, showing that the average distance between 
fragment centres is the same, irrespective of their size. The minimum distance 
(~ 8 fm for the smallest neighbour fragments) is the same in all cases. The 
points corresponding to big fragments {ZiZj > 700) are close to the hmit fine 
(calculated for two equal touching spheres at normal nuclear density) for the 
lighter system. It is not the case for the heavier system (middle panel). 

Moreover fig. 8 can also provide a rough estimate of the freeze-out volume. dZ 
is the infinitesimal number of charges in the volume element dT^ = Anr^dr. 
The upper panel in fig. 8 corresponds to Xe+Sn single source at 200 fm/c. 
The distribution has a nearly Gaussian shape, of mean f and full width at 
half maximum FWHM. A sphere of radius f + FWHM/2, indicated by the 
arrow in the figure, englobes most of the fragment centres and has a volume 
of ~ 2.8 Vq] the middle panel concerns the Gd+U single source at 240 fm/c: 
the volume of the sphere of radius f + FWHM/2 is ~ 7.9 Vq. Volumes at 
normal density Vq = (l-2)^Ato( fm^'^ are calculated for the masses of the 
sources given in table 1. The large difference between the volumes of the 
two systems, which may look surprising, is coherent with the infiuence of the 
system size on the Coulomb and radial expansion [50] , as already suggested by 
the average kinetic energies of the fragments (figs. 2 and 3). It might even be 
more pronounced, taking into account the underestimation of the compactness 
of the Xe+Sn system visible in the small discrepancy between the calculated 
and experimental correlation functions. The difference in the volume of the 
Xe+Sn source within 40 fm/c around the average freeze-out instant (figure 8 
and table 1) gives an estimate of the uncertainty on the determination of 
this volume. Note that the volume obtained at 240 fm/c for Xe+Sn compares 
well with that estimated by a simpler simulation using the measured charge 
partitions [54]. 



5 Conclusions 

Energy spectra of different fragments, issued from central collisions of 32 AMeV 
^^^Xe projectile on a "^"*Sn target and 36 AMeV ^^^Gd projectile on a "■"*U tar- 
get were presented. A detailed study of the evolution of the average energy 
versus fragment charge, and of the average sizes of the ordered fragments ver- 
sus fragment multiplicity evidences the particular role of the largest fragment 
in each event. Experimental relative velocity correlation function of different 
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types were built. Similarities between the two systems were inferred for the 
intra- fragment distances, while the position of the larger fragments relative 
to the centre of mass seemed to depend on the size of the multifragmenting 
source. 

All experimental results were satisfactorily well described by a stochastic 
mean-field model calculation, proving the capability of this dynamical ap- 
proach to reproduce in detail the experimental data and not only global, 
average observables. The observed differences noted for the Xe-|-Sn system 
may reveal that the topology of the events at freeze-out is not identical to 
the experimental one: fragmentation products seem to be more uniformly dis- 
tributed in space in the experimental configuration. Conversely the agreement 
between simulation and experiment is excellent for Gd+U, allowing to state 
that the simulation provides plausible reaction scenario and freeze-out topol- 
ogy. In the BOB image, involving the spinodal decomposition of the source 
(about 100 fm/c after the beginning of the collision), the freeze-out instant 
may be seen as the end of a time interval in which all the fragments are well 
separated, i.e. that moment where the last two closest fragments move away at 
~2 fm. This fact is on average accomphshed at 200 fm/c from the beginning 
of the simulations for the lighter system and at 240 fm/c for the heavier one. 
At the same moment: 240 fm/c, when the fragment multiplicity is definitively 
stable in both reactions, the radial fragment distribution refiects the infiuence 
of the system size on the Coulomb and radial expansion. The freeze-out vol- 
umes estimated at this moment are several times bigger than the volume Vq 
of the single source calculated at normal nuclear density. The higher compact- 
ness of the lighter system was quantitatively confirmed, as well as the more 
central position of the larger fragments. 
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